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ABSTRACT 



We report on the detailed modelling of the X-ray spectra of three likely neutron stars. The 
neutron stars, observed with XMM-Newton are found in three quiescent X-ray binaries in the 
globular clusters: w Cen, M 13 and NGC 2808. Whether they are accreting at very low rates or 
radiating energy from an accretion heated core, their X-ray spectra are expected to be those of 
a hydrogen atmosphere. We use and compare publicly available hydrogen atmosphere models, 
with constant and varying surface gravities to constrain the masses and radii of the neutron stars. 
Thanks to the high XMM-Newton throughput, and the accurate distances available for these 
clusters, using the latest science analysis software release and calibration of the XMM-Newton 
EPIC cameras, we derive the most stringent constraints on the masses and radii of the neutron 
stars obtained to date from these systems. A comparison of the models indicate that previously 
used hydrogen atmosphere models (assuming constant surface gravity) tend to underestimate 
the mass and overestimate the radius of neutron stars. Our data constrain the allowed equations 
of state to those which concern normal nucleonic matter and one possible strange quark matter 
model, thus constraining radii to be from 8 km and masses up to 2.4 Mq. 



Subject headings: Stars: neutron — Dense matter 
ij Cen, M 13, NGC 2808) — X-rays: stars 

1. Introduction 

Forty years after their discovery, the nature of 
the material making up the core of neutron stars 
remains largely unknown. At densities above a few 
times the equilibrium density of nuclear matter, 
models predict the existence of exotic components 
such as pio n or kaon condensates or unconfined 
quarks (e.g. iLattimer fc Prakashll2007l) . The ex- 
citing idea that neutron stars may contain exotic 
forms of matter makes them of prime interest, not 
only for astrophysics but for physics in general. 

Different equations of state of dense matter pre- 
dict different maximum masses and different mass- 
radius relationships. So far accurate mass mea- 
surements have been made for radio pulsars giv- 
ing masses ranging from 1.18±[] |]2 M0, for the 
case of one of the pulsars in PSR J1756-2251 
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([Faulkner et all 120051) to 1.4414±0.0002 M .^^. for 
PSR B1913-H6 IWeisberg fc Tavloill2005l) . Such 

low values can be accommodated by a wide range 
of equations of states and do not provide a strong 
constraint on the composition of dense matter. 
There is however growing evidence that neutron 
stars, as massive as 2 Mq, may also exist, in 
particular in accreting X-ray binaries, as inferred 
from the study of the properties of kilo-Hz QPOs 
(Quasi-Periodic Oscillations) (e.g. the neutron 
star in 4U 1636-53 6 which is estimated to have a 
mass of 1.9-2.1 M^ ISarret. Olive fc MIIIctI (|2005f) ) 
and binary millisecond pulsars that are no longer 
accreting (e.g. PSR J0751-I-1807, which has a 
mass measured through relativi stic orbital decay 
of 2.1±0.2 Mq^. lNice et"al]l2005[ ). If the latter es- 
timates are confirmed (by dynamical mass esti- 
mates of the binary components), they would rule 
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out some exotic forms of matter, such as quarks 
or pion condensates. 

While there are some accurate mass measure- 
ments aheady available, the situation is different 
with radii which are much harder to constrain. X- 
ray spectroscopy is one promising avenue to fol- 
low to tackle this issue. So far, there has been 
one single (and still debated) measurement of red- 
shifted absorption lines in the combined early and 
late type I X-ray burst spectra of the X-ray bi- 
nary EXO 0748-676. The line energies are con- 
sistent with FeXXVI Ha (n = 2-3) and FeXXV 
Heof (n = 2-3), respec tively, both at the same red- 
shift z = 0.350±0.005 dCottam . Paerels fc Meildei 
2OO2I ). The measured redshift provided a direct 



several neutron star systems, all providing plau- 
sible values fo r the neutron st ar radiu s (typically 



2002; Gendreetal 



estimate of Mass/Radius (M/R) {M/R = ^(1 - 
(l-|-z)~^). Assuming a canonical mass of 1.4 Mq, 
this would imply a radius of about 9 km. Fol- 
lowing this result, and using the burst properties, 
Ozell (12006*) estimated both mass and radius sep- 



arately. This led to a massive neutron star with a 
mass M > 2.10 ± 0.28Mo (and R> 13.8 ± 1.8 km) 
in EXO 0748-676, a result suggesting again that if 
this system is typical, exotic forms of matter, such 
as condensates and unconfined quarks do not exist 
in neutron star cores. 

Another promising way of inferring radii is 
through observations of quiescent X-ray emission 
from neutron stars for which the distance (d) 
can be estimated reliably [Foa — {Roo/dYaT^, 
where Foo, Too are the flux and radiation tem- 
peratures redshifted to the Earth and i?oo = 
R/ \Jl — 2GM /Rc^ is the radiation radius). Whether 
the energy reservoir is the heat deposited deep in 
the neutron st ar crust during the outburst phase o f 
the transient ( Brown. Bildsten fc Rutledge 1998f ). 
or sustained by a low- l evel of radial accretion 
( Van Paradiis et "al]ll987t ) (via an advection dom- 
inated accretion flow), the X-rays originate from 
the atmosphere of the neutron star. The spectrum 
radiated by the atmosphere will depend on its ac- 
tual composition, the strength and the structure 
of the magnetic field. In general, it is thought 
that the old neutron stars in globular clusters 
have low magnetic fields and hydrogen-rich atmo- 
spheres (gravitatio nal settling of heavier elements 
occurs rapidly, see lRutledge et al. I (|2002l) and ref- 
erences therein). Early non magnetic hydrogen 
atmosphere models have been shown to provide 
adequate fits to the quiescent X-ray spectra of 



around 10 km. lRutledge et al 
2003al lbh. Compared to neutron stars in the field. 



those in globular clusters have accurate distance 
estimates, and the fitting of their very soft spectra 
is made easier by the fact that at least the value of 
the interstellar absorption derived from the opti- 
cal extinction is well known, even if the absorption 
intrinsic to the system is largely unknown. 

Following the work done bv lHeinke et al. I (|2006h 
on Chandra data of the quiescent neutron star X- 
ray binary X7 in 47 Tucanae, we apply recently 
improved neutron star atmosphere models, avail- 
able in the latest release of th e XSPEC spec tral 
fitting package (version 12.3.0, Arnaud 1996f) . to 
three neutron stars in quiescent X-ray binaries 
that we have observed in three globular clusters 
with XMM-Newton. This is part of an on-going 
program to locate, identify and cl assify faint X-ray 
sources in glo bular clus ters (iServillat et al.l 2007 : 
Gendre et al] l2003al lbl: I Webb et all l2004 l2006^ . 



with an aim to determining the internal energy 
source that slows down the inevitable collapse 
of globular clusters. Early fitting of the X-ray 
spectr a with the first hydrogen-atmosphe re mod- 
els of iZavlin. Pavlov, fc Shibanov ( 19961 ) assum- 
ing a constant surface gravity (log gs=14.385, cor- 
responding to a 1.4 Mq and 10 km radius neu- 
tr qn star), were report ed for uj Cen and M 13 



in Gendre et al.l (2003al3). However, as empha- 
sised by Heinke et al.l (| 2006^, using models with 



appropriate surface gravity for each fitted value of 
the mass and radius of the neutron star is impor- 
tant when interpreting high quality X-ray spectra, 
which is indeed the case for the spectra measured 
with XMM-Newton. Applying the improved mod- 
els is the main motivation of this paper, which fur- 
ther benefits from improved data analysis software 
and calibration of the EPIC instruments compared 
to the earlier published results on these sources. 

In the next section, we describe the observa- 
tions, recalling the essential parameters of the 
clusters (distance, extinction,...), the data reduc- 
tion and present the results of the spectral fitting. 
We use, for comparison purposes three different 
models, described hereafter. 
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2. Observations and data reduction 

We have observations of the three hkely neu- 
tron stars in three different globular clusters, lo 
Centauri (cj Cen), M 13 and NGC 2808. Observa- 
tions of each of these clusters were made with the 
X-ray observatory XMM-Newton. All three EPIC 
cameras were used in the full-frame mode with the 
medium filter. Further information about these 
observations can be found in Tabled] 

To reduce these data we used the latest ver- 
sion of the XMM-Newton Science Analysis Soft- 
ware (SAS, version 7.0). This has many improve- 
ments over the earlier versions of the SAS used 
to reduce lo Cen and M 13 (version 5.3.3) such as 
upgraded EPIC calibration, resulting in a much 
better cross calibration among the EPIC instru- 
ments and which includes modelling of spatial and 
temporal response dependencies. Improvements 
to the bad pixel finding algorithms (decreasing the 
noise at low energies significantly, important for 
neutron stars which emit mainly at low energies), 
the vignetting correction and exposure maps have 
also been made (see the SAS release noteqll). The 
MOS data were reduced using the 'emchain'. The 
event lists were filtered, so that 0-12 of the prede- 
fined patterns (single, double, triple, and quadru- 
ple pixel events) were retained and the high back- 
ground periods were identified by defining a count 
rate threshold above the low background rate and 
the periods of soft proton flares were then flagged 
in the event list. We also filtered in energy. We 
used the energy range 0.2-10.0 keV, as recom- 
mended in the document 'EPIC Status o f Cal i- 

2002h . 



bration and Data Analysis' (jKirsch et al 



The pn data were reduced using the 'epchain' of 
the SAS. Again the event lists were filtered, so 
that 0-4 of the predefined patterns (single and dou- 
ble events) were retained, as these have the best 
energy calibration. We again filtered in energy, 
where we used the energy range 0.2-10.0 keV and 
we also filtered for the soft proton flares. A sum- 
mary of the observations and good time intervals 
for each observation and camera is given in Ta- 
ble ffl 

We extracted the w Cen spectra using circles of 
radii ~45" (to include at least 90% of the avail- 
able flux from the neutron star) centred on the 




Fig. 1. — Image showing the region used to extract 
the M 13 neutron star spectrum (brighter source) 
and the region excluded due to the source found 
at 18.4" from the neutron star, superposed on the 
MOS 1 data (0.2-10.0 keV). 



source. For the neutron star in M 13 we used an 
extraction radius of ^25" due its close proximity 
to other sources and excluded a region of radius 
15" around a second source that fell in the extrac- 
tion region (the large region contains at least 80% 
of the total source counts from the neutron star 
and less than 1 0% of the counts from the neigh- 
bouring source, Ehle et al. ( 2006f l). see Fig[T] Due 
to over-crowding in the centre of NGC 2808 we 
used extraction regions of 8", which includes 50% 
of the total sour ce counts from the neutron star 
(jEhle et allbood ). see FigH We used a similar 
neighbouring region, free from X-ray sources to 
extract a background spectrum. We rebinned the 
MOS data into 15 eV bins an d the yn d a ta int o 



5eV bins as recommended in Ehle et al 



(l2006l ). 

We used the SAS tasks 'rmfgen' and 'arfgen' to 
generate a 'redistribution matrix file' and an 'an- 
cillary response file', for each spectrum. 

3. Spectral Analysis 

We exploit three different hydrogen atmo- 
sphere models that are publicly available in Xspec 
12.3.0. These are the basic neutron star atmo- 
sphere model (nsa) that we used previously in 



Gendre et al 



(|2003allbl '). This model includes a 
uniform surface (effective) temperature, either a 
non-magnetised neutron star or with a field B = 
10^^ or B = 10^^ G and a radiative atmosphere 
in hydrostatic equilibrium. The nsagrav model 
used is similar, but allows for a variety of sur- 



i |http://xmm.vilspa.esa.es/sas/7.0.0/documentatioii/releasenotes/| face gravitational accelerations, ranging from 10 
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Table 1: Summary of the globular clusters (GC) and their observations. 





Date 


JJet. 


-L ohs 


i 1 


UlSt. 




u) Cen 


zUUi Aug 13 


iVll 


40 


40 


5.3 


U.Ob/ 






M2 


40 


40 










pn 


40 


40 






M 13 


2002 Jan 28 


Ml 


17 


12 


7.7 


0.011 






M2 


17 


11 










pn 


14 


8 






M 13 


2002 Jan 30 


Ml 


18 


14 










M2 


18 


14 










pn 


13 


8 






M 13 


1990 Jun 1 


PSPC 


46 


46 






NGC 2808 


2005 Feb 2 


Ml 


41 


41 


9.6 


0.128 






M2 


41 


41 










pn 


40 


31 







Dct. indicates the detector used (Ml/M2= MOS 1 or 2), Tohs gives the observation time in ks and GT I gives the good timing 
interval (in ks) after filtering for soft proton flares. The distances (Dist.) arc taken from Harris (1999) and are given in kpc. 
The column densities (n^, xlO^^ cm~-^) to these clusters are from iHarrisi (I199& ): iPredehl fc Schmitli (il995i ). 




of this source in M 13 (jVerbunll l200ll ) as m 
Gendre et all (|2003bf). as the ROSAT data ex- 
tends below the XMM-Newton data, down to 0.1 
keV (only 11.5 ks of data exists in the archives for 
uj Cen, which is insufficient to improve our spectra. 
No data exist for NGC 2808). The ROSAT PSPC 
has a poorer angular resolution than the EPIC 
cameras (PSPC = 25" and MOS = 6") and the 
spectrum extracted contains data from both the 



Fig. 2.— MOS 1 image (0.2-10.0 keV) showing 
the region used to extract the NGC 2808 neutron 
star spectrum and it's proximity to the bright cat- 
aclysmic variable. 

to 10^^ cm s~^, adapted to the masses and radii 
investigated. Finally we use th e nsatmos model 



described bv lHeinke et al. ( 20061 ). This model in 



eludes a range of surface gravities and effective 
temperatures, and incorporates thermal electron 
conduction and self-irradiation by photons from 
the compact object. It also assumes negligible 
(less than 10^ G) magnetic fields and a pure hy- 
drogen atmosphere. 

For each neutron star, we fitted the spectrum 
obtained from all three cameras simultaneously. 
However, for the neutron star in M 13 we also in- 
cluded the ROSAT PSPC archival observations 



the neutron star, see Fig. [TJ We have extracted 
the XMM-Newton EPIC MOS spectrum of this 
source. The best fit to this spectrum is either 
a power law, r=1.75±°;^J, Cstatistic = 5.96 (8 
bins) or a bremsstrahlung model, kT=4.49±2"g(,'''', 
Cstatistic = 5.98 (8 bins). This source contributes 
approximately 28% of the counts in the ROSAT 
band (0.1-2.4 keV). To take this source into ac- 
count, we have allowed the normalisation for the 
ROSAT data to float, to account for the uncertain 
cross-calibration between the two observatories. 

All the data were binned to contain a mini- 
mum of 20 counts per bin where possible and we 
used the technique to judge the goodness of 
the model fits. However, for lo Cen and NGC 
2808 the data were binned with 15 counts per 
bin (to increase the number of data points). In 
these cases we used the C-statistic to judge the 
goodness of the fit, as the number of counts per 
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Table 2: Results of the spectral fitting for the three neutron stars with the three different hydrogen 



atmosphere 


models, as described in 


Section [3l 








Cluster 


NSA 


NSAGRAV 


NSATMOS 




:> ° 1 


w Cen 


;o5(T)=5.99±-™ 
M=1.76±'i:lZ 

R=n.30±:f„ 

6. 3 Op 

Cstat= 91.28, 85 dof 


7V^=0.11±°:- 
log{T)=5.99±ll° 

i?=10.00±I;- 
Cstat= 91.26, 85 dof 


/o5(T)=5.98±-- 
M=1.66±-- 
i?=11.66±I;°^ 
Cstat= 91.35, 85 dof 


4.91 X 




M 13 


Zo5(T)-6.00±°;- 
M=1.38±°«« 
i?=9.95±°^* 

Xf, = 1.10, 62 dof 


iV^=0.013±-°- 
logiT)=6m±l°t 
Af=1.39±°-^^ 
i?=9.95±^^^ 
X', = 1.10, 62 dof 


log{T)=6.00±rol 
M=1.30±°;« 
i?=9.77±«;- 

X' = 1.08, 62 dof 


5.08 X 




NGC 2808 


7Vff=0.17±«;- 
log{T)=6M±lZ 
M=0.67±°ll 

Cstat=17.59, 19 dof 


;o5(r)=6.08±°-5'« 

0.65p 

R=7A8±lil 
Cstat=18.12, 19 dof 


NH^0A6±'>lt 
log{T)=6m±lll 

0. 41p 

i?=6.10±";^^ 

1 . lOp 

Cstat=16.95, 19 dof 


1.02 X 





The column density, Njj is xlO^^ cm~^, the temperature (T) is given as the logarithm of the temperature in Kelvin, the masses 
are in solar units and the radii in kilometres. All errors are given at the 90% confidence limit for the one interesting parameter. 
To calculate the errors, the mass and distance were held fixed and all other parameters were allowed to vary, except when 
calculating the error on the mass, when the radius was held fixed. A 'p' after the error value indicates that the hard limit of 
the model was reached. The estimated unabsorbed bolometric luminosity for each neutron star is also given. 



bin is no longer strictly in the Gaussian statis- 
tics regime but approaches that of the Poisso- 
nian statistics regime. Even though we are at the 
limit between the two regimes, the C-statistic has 
been shown to work well at even higher counts 
(|Nousek fc Shudll989l ) and should therefore give 
good results. T h e inve rse is not necessarily true. 
Nousek fc Sh"u3 ()l989l) showed that using the 



Levenberg-Marquardt algorithm (the algorithm 
used in Xspec) for small numbers of events per 
bin gives a systematic bias. 

3.1. Absorption in the spectra 

For each model we included (photoelectric) ab- 
sorption along the line of sight to the object. This 
was initially fixed at the value determined for each 
cluster (see Table [Ij and then allowed to vary. In 
every case we found that we required additional 
absorbing material, assumed to be additional gas 
intrinsic to the system. This was typically 30- 
60% more i.e. w Cen required an n^f-^^O.ll xlO^^ 
cm~^, an increase of 60% with a ftest probability 
of 0.0099, indicating that it is reasonable to in- 
clude additional material. The normalisation was 



fixed to the value corresponding to the distance to 
the cluster (see Tabled]) and the masses and radii 
were initially frozen to 1.4 Mq and 10.0 km, but 
then allowed to vary. 

We then tried to constrain the nature of the 
additional absorbing material, in the same way 
as Heinke et al. (2006), by employing the Xsyec 
model vphabs and using the abundances corre- 
sponding to those of each of the clusters. For M 
13 we chose the iron abundance to be 3% solar 
([Fe/H] = -1.5), the abundances of C, N and O to 
be 5% solar ([X/H] = -1.32), the abundances of Ne 
to Ca to be 4% solar ([X/H] = -1.38) and that of 
helium to be of solar abundance, following Cohen 
& Melendez (2005). For w Cen the situation is 
more complicated as there are three distinct pop- 
ulations of stars in this cluster, see e.g. Origlia 
et al. (2003). We elected the abundances from 
the largest population in the cluster, basing our 
choice on the fact that the highest probability was 
that the donor star comes from the largest popula- 
tion. We adopt the iron abundance to be 3% solar 
([Fe/H] = -1.58), the abundances of C, N and O 
to be 6% solar ([X/H] = -1.24), the abundances of 
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Fig. 3. — Unfolded spectrum of the neutron star 
in M 13 with the best nsatmos model fit and the 
residuals to the fit. The different symbols indicate 
the data from the different instruments used: open 
diamonds show PSPC data, filled circles and tri- 
angles show MOS 1 data, filled squares and empty 
circles show MOS 2 data and filled diamonds and 
empty squares show pn data. 

Ne to Ca to be 6% solar ([X/H] = -1.19) and again 
that of helium to be of solar abundance, following 
Origlia et al. (2003) and Norris (2004). For NGC 
2808 we select the iron abundance to be 9% solar 
([Fe/H] = -1.06), the abundances of C, N and O 
to be 32% solar ([X/H] = -0.5), the abundances 
of Ne to Ca to be 16% solar ([X/H] = -0.8) and 
that of helium to be of solar abundance, follow- 
ing Castellani et al. (2006) and Gratton (1982). 
Again, however, it is believed that NGC 2808 has 
as many as three populations of stars, with dif- 
ferent abundances (Piotto et al. 2007). We chose 
the abundances for the only population with suffi- 
cient information in the literature. We then mod- 
elled the neutron star spectra fixing the photo- 
electric absorption (phabs) to be that of the clus- 
ter and allowing the photoelectric absorption with 
variable abundances (vphabs) to vary (using the 
abundances given above) . The results that we ob- 
tained gave very small masses and/or radii, often 
small enough to hit the lower mass or radius limit 
allowed in the models (see below for values). 



Fig. 4. — Unfolded spectrum of the neutron star 
in iv Gen with the best nsatmos model fit and the 
residuals to the fit. The different symbols indicate 
the data from the different instruments used: filled 
circles show MOS 1 data, filled squares show MOS 
2 data and filled diamonds show pn data. 



0.003 
0.002 
0.001 
0.000 
-0.001 
-0.002 



Energy (keV) 

Fig. 5. — Unfolded spectrum of the neutron star 
in NGC 2808 with the best nsatmos model fit and 
the residuals to the fit. The different symbols used 
are the same as those employed in Fig. |4l 

We thus decided to determine whether our data 
are of sufficiently good quality to constrain the na- 
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ture of the additional absorbing material. We sim- 
ulated a spectrum similar to that obtained for lo 
Cen, using the simple phabs model. We modelled 
this spectrum in Xspec using first the phabs model 
and then the vphabs model. We obtained similar 
Xl values (1.17, 78 dof when using phabs and 1.29, 
78 dof when using vphabs), but again the masses 
and radii were very low when the vphabs model 
was employed, only ~80% of the value given to 
simulate the spectrum as opposed to values within 
10% when using the phabs model. We also sim- 
ulated the same spectrum but using the vphabs 
model. Again, modelling this spectrum in Xspec 
using the vphabs model we found values as low 
as 33% of the original value, yet when modelling 
with the phabs model our values are within 5% 
of our original values. Further, we obtain equally 
good fits using vphabs to our simulated spectrum 
(x^=1.09, 78 dof) when we choose abundances 
that are different by 300% to those used in the 
simulated spectrum, and we still have the prob- 
lem that either the mass or the radius is very small. 
We therefore conclude that the quality of our data 
is insufficient to constrain the nature of the addi- 
tional absorbing material and we use the phabs 
model only. 

3.2. Contribution from neighbouring sources 

We also investigated whether we required an 
extra power law tail, especially for the two neutron 
stars found in proximity of other sources i.e. the 
neutron stars in M 13 and NGC 2808. This is 
important as substantial flirx from other sources 
is likely to be the dominant source of counts in 
the higher energy bands (i.e. above 1 keV). This 
may tend to systematically bias the spectral fits 
to higher temperatures, and thus smaller radii. 

To do this we extracted the spectra of the 
neighbouring sources and determined their pho- 
ton indices (r-1.75 for M 13 and r-1.5 for NGC 
2808). We then estimated the contribution from 
the neighbouring sources in the extraction region 
for each neutron star (8% for M 13 and 20% for 
NGC 2808, see Figs [T] and [2]). We used these 
values and refitted our spectra. Interestingly, we 
found values very similar to those obtained when 
no power law was included. Even allowing the 
normalisation to increase, increasing the weight of 
the power law, made little difference to the fits. 
Ftests show that adding the power law is not nec- 



essary. Probability results for fitting nsa, nsatmos 
and nsagrav models with and without the addi- 
tional power laws to the neutron star in M 13 are 
0.31, 0.22 and 0.199 respectively and to the neu- 
tron star in NGC 2808 are 0.13, 0.06 and 0.08 
respectively. We therefore conclude that adding a 
power law is not required by the data. 

The spectra of the three neutron stars can be 
found in Figs.[3][5l along with the best fitting nsat- 
mos model. 

3.3. Modelling the spectra 

The principal goal of our spectral fitting is 
to self-consistently constrain the allowed space in 
mass and radius using our three neutron stars, 
but also to test the reliability and accuracy of 
the three models examined. Table [2] gives the 
best fits for each neutron star fitted with each 
of the three models, along with the best fitting 
values for the column densities, surface tempera- 
tures of the neutron stars, masses and radii. Er- 
rors are also given at the 90% confidence limit for 
the one interesting parameter. To calculate the 
errors, the mass and distance were held fixed and 
all other parameters were allowed to vary, except 
when calculating the error on the mass, when the 
radius was held fixed instead of the mass. We 
also give the estimated bolometric luminosity of 
the neutron star. We found that we did not re- 
quire any additional parameters, such as lines or 
edges to fit these data. However, this may be 
due to the quality of the data, where deeper ob- 
servations may reveal spectral features that will 
be useful to constra in the gravitational redshift 
at the NS surface (IBrown. Bildsten fc Rutledge 



19981: iRutledge et al.ll2002l: iHeinke et al.ll2003l ) 

We use the steppar command in XSPEC to vary 
both the radius and mass parameters simultane- 
ously, allowing the temperature to vary as well to 
find the best fit. The minimum (and maximum) 
radii allowed with these models are 5 km (30/20 
km) for the nsatmos and nsa models respectively 
and 6 km (20 km) for the nsagrav model. We 
chose an inferior limit of 8 km and a superior limit 
of 18 km as all three models gave stable results 
when fitting the data in these regions. The min- 
imum (and maximum) masses allowed with these 
models are O.SM© (3.O/2.5M0) for the nsatmos 
and nsa models respectively and O.3M0 (2.5Mq) 
for the nsagrav model. We chose the region be- 
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Fig. 6. — Contour plot showing the results of 
modelling the neutron star in M 13 with the 
xspec models: nsa (solid/red line), nsatmos (dot- 
ted/blue line) and the nsagrav (dashed/green 
line). The 90% and 99% confidence contours for 
the neutron star masses (M©) and the radii (km) 
are plotted along with neutro n star e quations of 
state taken from .Lattimer fc PrakashI (12007,) 



Fig. 7. — Contour plot showing the results of mod- 
elling the neutron stars in lo Cen (solid/red line) 
and NGC 2808 (dotted/blue hue) with the xspec 
model nsatmos. The 90% and 99% confidence 
contours for the neutron star masses (M©) and 
the radii (km) are plotted alo ng with neutron star 
equations of state taken from lLattimer fc Prakash 
^2001) 



tween 0.5 and 2.3M0 as again all three models 
gave stable results in these regions. We show 90% 
confidence, and 99% confidence (Ax^ = 4.61, and 
9.21 respectively) contours in neutron star mass 
and radius. Fig. [6] shows the results of modelling 
the neutron star in M 13 with all three neutron 
star atmosphere models. This figure indicates 
that the two models, nsatmos and nsagrav, give 
comparable results, whereas the nsa model gives 
differing results. The nsa model is more accu- 
rate in constraining Rqo than the nsatmos and 
nsagrav models, due to a degeneracy in spectral 
shape variations between the surface gravity and 
the surface temperature. However, here we calcu- 
late a range of neutron star masses and radii, for 
which the nsatmos and nsagrav, thanks to their 
variable surface gr avities, are better a dapted (see 



the discussion in iHeinke et all l2006f) . The use 



of these fixed gravity models for testing neutron 
stars with a variety of masses and radii (and there- 
fore gravities) is not therefore strictly appropriate. 



This is borne out by the modelling. For this rea- 
son and to simplify the following plots, we ignore 
the nsa model (as it is less accurate) and we plot 
only the nsatmos model (as the nsagrav and nsat- 
mos models are similar) to show the constraints 
on mass and radius determined by modelling uj 
Cen and NGC 2808 (Fig. [T]). The loci of mod- 
els for the equati ons of state for dense matte r are 
those described in Lattimer fc PrakashI (|2007l ) and 
Lattimer fc PrakashI ( 200ll) which include diverse 



Prakash. Cooke fc Lattimer 



equat ions such as: SQM ■ 
(Il995l ). a Strange Quark Matte r mod el, 



PAL 



Prakash. Ainsworth fc Lattimer (|l988h a neu 



tron and pro ton model using a schemati c poten- 



tial, CM - iGlendenning fc Moszkowskil (|l99ll ). 



a model containing neutrons, protons and hy- 
perons using a field theoretical appro ach, GS - 
Glendenning fc Schaffner-Bielich a model 



containing neutrons, protons and kaons using a 
field theoretical approach. 
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Fig. 8. — Contour plot showing the results of mod- 
elling the neutron stars in M 13 (dotted/blue line), 
uj Cen (solid/red lin e) and X7 in 47 Tu c (dashed- 
dotted/green line) (jHeinke et al.l l2006l ) with the 
xspec model nsatmos. The 99% confidence con- 
tours for the neutron star masses (Mq) and the 
radii (km) are plotted along with neutron star 
equat ions of state taken from .Lattimer fc PrakashI 
The grey hashed region indicates the re- 
gion not investigated with the models. 



4. Discussion and Conclusions 

As seen in Section |3] even if the nsa model 
is more accurate in constraining Roo than the 
nsatmos and nsagrav models, when calculating a 
range of neutron star masses and radii, we find 
that the nsatmos and nsagrav models are better 
adapted. We conclude that the use of fixed grav- 
ity models for testing neutron stars with a vari- 
ety of masses and radii (and therefore gravities) is 
not strictly appropria te as previously indicated by 



not strictly appropria te as previously mdicatea i 
Gansicke etall (|2002h and lHeinke et al.l (|2006l ). 



Using the res ults from Fig s . El an d [71 and the 
results found bv lHeinke et al. (|2006l ) when fitting 
Chandra observations of X7 in 47 Tuc, we show 
the allowed equations of state in Fig. |8l Modelling 
the neutron star in M 13 alone shows that the data 
do not favour the stiffer equations of state, such 
as MSO-2 and PAL 1. We combine these results 



with those from modelling the neutron star in lo 
Cen. The low quality data for the neutron star in 
NGC 2808 results in very poor constraints for this 
object and hence we disregard this source in the 
discussion, although the results seem to be similar 
to those for the neutron star in M 13. The equa- 
tions of state that are satisfied by all three neutron 
stars fall in the middle of the diagram and includes 
the equations of state of normal nuclconic matter 
and one strange quark matter model. The equa- 
tions allowed are GSl, PAL 6, AP3 & 4, GM3, 
FSU, SQM3, ENG and MPAl, with radn above 8 
km and masses up to 2.4 M©. 

This article was based on observations obtained 
with XMM-Newton, an ESA science mission with 
instruments and contributions directly funded by 
ESA Member States and NASA. The authors 
thank J. Lattimer for providing them with the 
mass-radius relationships for the different equa- 
tions of state shown in Figures [MHl and B. Gendre 
for help with the ROSAT spectrum. The authors 
also acknowledge the CNES for its support in this 
research and thank the (anonymous) referee for 
very valuable comments that have enabled us to 
improve the quality of this paper. 
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